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How to extract physics from HBT radius parameters 
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I review recent progress in the understanding of the connection between the space-time structure 
of the particle emitting source and the form of the two-particle correlation function in momentum 
space. Based on a new scheme for calculating the HBT radius parameters from the emission 
function, strategies are suggested to separate for rapidly expanding sources the information 
on the spatial and temporal structure of the source. To this end a new fitting function for 
the two-particle correlation function is proposed. Its usefulness is demonstrated for a typical 
expanding model source, and it is shown how the dependence of the resulting fit parameters 
on the momentum of the particle pair can be used to measure the longitudinal and transverse 
expansion of the source. 

1. INTRODUCTION 

2 . 

The only known way to obtain direct experimental information on the space-time struc- 
ture of the particle emitting source created in a relativistic nuclear collision is through 
two-particle intensity (Hanbury-Brown-Twiss (HBT)) interferometry [ |fl. The goal of 
this method is to extract the space-time structure of the source from momentum spectra 
which are the only measurable quantities, making use of the quantum statistical correla- 
tions between pairs of identical particles. This information is crucial for an assessment of 
theoretical models which try to extract the energy density of the source from the measured 
single particle spectra and particle multiplicity densities in momentum space. Reliable 
estimates of the source volume and the energy density are, on the other hand, indispens- 
able for an experimental proof that high energy collisions can successfully generate large 
volumes of matter with extreme energy density, where a transition into deconfined quark 
matter might be possible. 

For many years HBT interferometry of hadron-hadron and nucleus-nucleus collisions 
was motivated by the naive expectation, based on the experience with photon interfer- 
ometry for stars, that the width of the two-particle correlation function in the relative 
momentum directly measures the geometric size of the source. This expectation is wrong. 
Unlike stars in the universe, the sources created in hadronic or heavy- ion collisions may 
feature inhomogeneous temperature profiles and strong collective dynamical expansion. 
We now know [ |, |, H |, || that for such sources the HBT radius parameters generally 
don't measure the full source size, but only so-called "space-time regions of homogeneity" 
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[ 0, f|] inside which the momentum distribution varies sufficiently little so that the parti- 
cles can actually show the quantum statistical correlations. The size of these homogeneity 
regions varies with the momentum of the emitted particles, causing a dependence of the 
HBT parameters on the pair momentum [ g |, [§ |, |, 0, §, |, [10], [11], |g |13|, |TjJ. The 
detailed momentum dependence is, however, model-dependent, and in general it is not 
simple [0. The extraction of the strength of the collective flow from HBT data is further 
complicated by significant resonance decay contributions which also induce a momentum 
dependence of the HBT radius parameters and of the so-called "incoherence parameter" 

MM- 

The finite lifetime of the sources created in nuclear collisions leads to another compli- 
cation: the HBT radius parameters generally mix the spatial and temporal aspects of 
the source extension in a non-trivial and reference frame dependent way, in particular if 
at the same time the source undergoes collective expansion. The origin and pattern of 
this mixing was clarified in a recent series of publications from the Regensburg group [ 

O, [T2|, [T3L rj|. Two major new results have resulted from this work: (i) the discovery 



[ U of a new cross term in the two-particle correlation function mixing the outward and 
longitudinal components of the relative momentum between the two particles, and (ii) 
the new YKP (Yano-Koonin-Podgoretskii) parametrisation of the correlation function [ 



TT| , [13], [I4j . The latter permits the experimental determination of the longitudinal veloc- 
ity of the source volume element where most of the particle pairs originate (as a function 
of the pair momentum) and achieves a nearly perfect factorization of the longitudinal, 
transverse and temporal homogeneity regions of the source (again as functions of the pair 
momentum) in the source rest frame. Furthermore it provides for a clean separation of 
the longitudinal and transverse dynamics of the source. 

In this talk I review these new theoretical developments and exemplify them for a class 
of simple model emission functions for thermalized sources with collective transverse and 
longitudinal expansion and finite space-time geometry. 

2. EMISSION FUNCTION AND PARTICLE SPECTRA 

The single and two particle spectra are defined as 

Pi(p) = = E p (a+a p ) , (1) 

dN 

^2(Pi,P 2 ) = E X E 2 ^^ = E l E 2 {a+ l a+ 2 a P2 a pl ) (2) 

in terms of creation and destruction operators for on-shell particles with momenta pi, 
where (...) denotes an average over the source ensemble. Pi and P 2 are normalized to 
the average number of particles (N) and pairs (N(N — 1)) per event, respectively. The 
two-particle correlation function is defined as 

rf n n \- W *Mpi,p a ) 

C ^™~ {N{N-l))P^) Pl [&Y (3) 



For uncorrelated emission and in the absence of final state interactions [ [L6|] one can prove 
[ [L7|, [ll| a generalized Wick theorem for the factorisation of the 2-particle spectrum @) 
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and obtains 
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c( q ,K) = i± „:/\r', r (4) 



where q = Pi — P2 and K = (pi + P2)/2 denote the relative and total momentum of the 
particle pair, and the positive (negative) sign applies for bosons (fermions). Note that 
the second term is positive definite. 

These expressions can be further simplified and turned into a practical starting point 
for computations by introducing the emission function S(x,K). It is defined in terms 



of the classical source amplitude J(x) for creating a free pion state [ |17 via the Wigner 
transform of its associated density matrix 

S{X > K) = l 4k C ~ K ' V ^ ( X + ") J ( X " ") ) (5) 

and is the quantum mechanical analogue of the classical phase-space density which gives 
the probability for creating a free particle with four-momentum K at space-time point x. 
In terms of this emission function the single particle spectrum is given by 

Ek§^ = J d 4 xS(x,K) (6) 



where the r.h.s. is to be evaluated on-shell, i.e. at K° = E K = \/m 2 + K 2 . The two- 
particle correlation function is obtained from [ [19], |17|, 0, [H| 



where the r.h.s. must be evaluated at q = p% — P2, K = (p± + p2)/2 with pi on-shell. (This 
implies K-q = 0.) The approximation consists of replacing the single particle spectra at p\ 
and P2 in the denominator by the spectrum at K = (p± +p2)/2; it is exact for exponential 
momentum spectra and a good approximation in practice [ ||. The second equality in 
(0) defines a (independent) average over the emission function of which we will make 
abundant use below. A useful feature is that in (^) the emission function can be to very 
good approximation [ [| evaluated at K° = Ek, i.e. on the classical energy shell, since 
the typical source radii are larger than the Compton wavelengths of the observed hadrons. 
This warrants the replacement of the Wigner density S(x, K) by a classical phase-space 
distribution in practical calculations. 

Due to the on-shell constraint q-K = the four components of q are not independent, 
but related by 

q° = (3-q with /3 = -^A. (8) 

The Fourier transform in (|7|) is therefore not invertible, and the reconstruction of the 
space-time structure of the source from HBT measurements will thus always require addi- 
tional model assumptions. Furthermore, inserting (|]) into (0), iq-x = iq- (x — f3t), we see 
that the correlator C(q, K) mixes the spatial and temporal information in a non-trivial 
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way which depends on the pair velocity (3. Only for time-independent sources things 
become simple: the correlator then measures the Fourier transform of the spatial source 
distribution, however only in the directions perpendicular to (3 since the time integration 
leads to a 5-function S(/3-q). 

From Eq. (^) it is clear that, unless the emission function factorizes in x and K, 
S(x, K) = F(x)G(K) (in which case G(K) cancels between numerator and denominator), 
the correlator is a function of both q and K. If one parametrises it by a Gaussian in q (see 
Sec. |j) this results in K-dependent width parameters ("HBT radii"). In thermal sources 
x — K correlations which spoil such a factorization can be induced by temperature gradi- 
ents and/or collective expansion (with a 4- velocity u^{x))\ in both cases the momentum 
spectrum ~ exp[— p-u(x)/T(x)] of the emitted particles depends on the emission point. 

3. MODEL-INDEPENDENT EXPRESSIONS FOR THE HBT RADII 

One of the crucial questions is, of course, to what extent a measured K-dependence of 
the HBT radii allows for a quantitative reconstruction of the collective source dynamics. 
To answer it we must first learn more about the physical meaning of these "radii". To 



this end it is useful to write the emission function in the following form [ |TT|, [12|, [13 

1 



+ SS(x,K), (9) 



S(x, K) = N(K) S(x(K), K) exp 
where (with expectation values defined as in (0)) 

x M (K) = (x^), x M (K) = x tl - x M (K), (B^^K) = (x^). (10) 

This construction ensures that the term 5S has vanishing zeroth, first and second or- 
der moments and thus contains only higher order information on sharp edges, wiggles, 
secondary peaks, etc. in the source. It was shown numerically [ |12| to have negligible 
influence on the half width of the correlation function and to contribute only weak, es- 
sentially unmeasurable structures in C(K,q) at large values of q. Neglecting 5S, the 
two-particle correlation function (^) can be calculated analytically: 

C(K,q) = l + exp[-g^(5^)(K)] . (11) 

Please note that the point x M (K) of maximum emissivity at momentum K is unmeasurable 
[ O, Only the K-dependent effective widths ("lengths of homogeneity") (x^x u )(K) 
of the source of particles with momentum K are accessible by HBT interferometry. 
Actually, due to the on-shell constraint (||), only 6 linear combinations of the 10 vari- 



ances (a; A1 x l/ )(K) are measurable [ II]; in the case of azimuthal symmetry of the source 
around the beam axis, this number reduces to 4 out of 7. Which linear combinations oc- 
cur in practice depends on the way the correlation function is parametrised. The general 
form (|TT|) together with (H) still provide some freedom as to which components of q to 
keep as independent variables (see Sec. But whichever choice one makes, all the In- 
dependent parameters ( "HBT radii" ) in the resulting Gaussian function of q can be easily 
calculated from the variances (x^x"), i.e. by simple quadrature formulae, for arbitrary 
emission functions S(x, K). The relation between the HBT parameters and the variances 
is model-independent, i.e. it does not depend on the form of the emission function S(x, K). 



5 



4. STANDARD AND YKP FITS TO THE CORRELATION FUNCTION 



For the following discussion we employ the conventional [ [7|, || Cartesian coordinate 
system with z along the beam axis and K lying in the x-z-plane. The z-component of 
a 3- vector is labelled by / (for longitudinal), the x-component by o (for outward) and 
the ^/-component by s (for sideward). Then (3 S = such that q° = j3±q + fiiqi, with 
/3±_ = |KjJ/.fr being (approximately) the velocity of the particle pair transverse to the 
beam direction while fti is its longitudinal component. 

The standard Cartesian parametrization [ [| of the correlation function is obtained by 
using this condition to eliminate q° from Eq. (|TT|) . One obtains 



C(K, q) = 1 + exp 



- £ i£(K) ft <& 



i,j=s,o,l 



where the 6 HBT "radius parameters" Rij are given as [ |6|, pT 
4(K) = ((xi - - fot)) , i,j 



8,0,1, 



:i2) 



(13) 



through through the space-time variances of the source. For an azimuthally symmetric 
sample of collision events, C(q, K) is symmetric with respect to q s — ► — q s [ |TTJ. Then 
R 2 =]?2 o and 



C(K,q) 
i? s 2 (K) 



^(K) 



^(K) 



1 + exp [-i?2 (K )g s 2 - R 2 (K)q 2 - Rf(K)q? - 2R 2 ol (K)q oqi 

if), 

((x-{3j) 2 ), 

((5 -m 2 ), 

((5 - P±t)(z - A*)) • 



with 



(14) 

(15) 
(16) 
(17) 
(18) 



The cross-term (|18| ) was only recently discovered [ gf. Clearly these HBT radius param- 
eters mix spatial and temporal information on the source in a non-trivial way. Their 
interpretation in various reference systems, in particular the meaning of the generally 
non- vanishing cross-term R 2 ol , was extensively discussed in Refs. [ [51 O, [T3|, by analysing 
these expressions analytically for a large class of (azimuthally symmetric) model source 
functions and comparing with the numerically calculated correlation function (0). An 
important observation resulting from these studies is that the difference 



p2 



Rl-Rl 



(19) 

14] and thus provides access to the 
(t 2 ) — (t) 2 of the source [ ^ (more exactly: the duration of the particle 
emission process) . In heavy-ion collisions, due to rapid expansion of the source, one would 
generally not expect (t 2 ) to be much larger than either (x 2 ) or (y 2 ) (see however [ p3 



is generally dominated by the first term on the r.h.s. 
lifetime At = 



for possible exceptions near a phase transition to QGP). In the standard fit one is not 
sensitive to small values of At since Eq. ([yD then involves a small difference of two large 
numbers, each associated with standard experimental errors. The factor 0\_ < 1 in front 
of (t 2 ) further complicates its extraction, in particular at low K± where Ai(K) is usually 
largest (see below). 
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This problem is avoided in the Yano-Koonin-Podgoretskii parametrisation [ |24], |25], |TT 



13, |I~4|j of the correlation function for azimuthally symmetric systems. It is based on an 



elimination of q Q and q s in terms of q\ = q 2 Q + q 2 , q°, and q 3 in (|TT 



C(q, K) = 1 + exp [-R\ q\ - R 2 (qf - (q ) 2 ) - (R 2 + R 2 ) (q-U) 2 } , (20) 

with four K-dependent parameters R±, R\\, Rq, and If 1 where the latter is a 4-velocity 
with only a longitudinal spatial component: 

[/(K) = 7 (K) (1,0,0, <u(K)), with 7 = (1 - v 2 y 1/2 . (21) 

This parametrisation has the advantage that the fitted YKP parameters i?j_(K), R 2 (K), 
and -Rq(K) do not depend on the longitudinal velocity of the observer system. They (as 
well as v (K)) can be calculated from the variances (x^x u ) in any reference frame (see [ |T3| 
for explicit expressions), but their physical interpretation is easiest in terms of coordinates 
measured in the frame where v(K) vanishes. There they are given by [ |TT] 

Rl(K) = R 2 (K) = (f), (22) 
R 2 (K) = ((z-A5/^) 2 )-A 2 (y 2 )//5l^(5 2 ), (23) 

Rim = ((t-^ ± ) 2 )-(y 2 )//3i^(P), (24) 
where in the last two expressions the approximation consists of dropping generically small 



11 1 terms (for a quantitative discussion see [ ll|). The first expression (p2|) remains true 
in any longitudinally boosted frame. 

Eq. (0) shows that the YKP parameter i?o(K) measures directly (up to the neglected 
small terms) the time duration At(K) during which particles of momentum K are emitted, 
in the frame were the YKP velocity v (K) = 0. The advantage compared to the standard 
Cartesian fit is that here it is fitted directly, and no problems of differences of large 
numbers occur in its extraction. 



Since the standard Cartesian and YKP parametrizations (0) and ( p0|) of the correlator 
differ only by the choice of independent components of q, the two sets of HBT parameters 
must be related. One finds [ [U| 

R 2 S = Rl, (25) 
R 2 diS = R 2 -R 2 S = P 2 ±1 2 (R% + v 2 R\) , (26) 

R 2 = (l-p 2 )Rl + 1 2 ((3 l -v) 2 (R 2 + R 2 ) , (27) 
R\ = P± (-PiRj + 7 2 (A - vf {Rl + R\)) . (28) 

These relations provide a powerful consistency check on the experimental fitting procedure 
of the correlation function, of similar value as the relation [ [U], |1^] lim^ ± ^o(-Ro(K) — 
R S (K)) = which results from azimuthal symmetry. 
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5. A SIMPLE SOURCE MODEL 



For a quantitative discussion of the physical behaviour of the HBT radius parameters, 
in particular of their K-dependence, we use a simple model for the emission function of a 
finite expanding thermalized source [11]: 



S(x,K) 



M± cosh(r/ - Y) 



(27r)V2vr(Ar) 



exp 



K ■ u(x) 
T 



T 



TO 



r rj — t)q 
2(Ar) 2 ~ 2R 2 " 2(A V y 



(29) 



Here r = y/x 2 + y 2 , the spacetime rapidity rj = §ln[(t + z)/(t — z)] and the longitudi- 



parametrize the spacetime coordinates x^, with measure 



nal proper time r 
d A x = r dr drj r dr c 
longitudinal and transverse components of the pair momentum K. 



F = |ln[(l + A)/(l-A)]andM i 



+ K\ parametrise the 



E 

*j 3 
cr 
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Fig.l. The standard Cartesian parameters R s (a), R Q (b), Ri (c), and 
R^ (d) in the CMS for pion pairs with cm. rapidity Y = 1.5, as functions 
of Mi for 3 different values for the transverse flow rjf. The thick lines 
are exact numerical results from Eqs. ( |15| - [18|) , the thin lines are obtained 
from the analytical approximations given in Ref. [|5[]. (Figure taken from 
Ref. 



T is the freeze-out temperature, R is the transverse geometric (Gaussian) radius of the 
source, To its average freeze-out proper time, At the mean proper time duration of par- 
ticle emission, and Arj parametrises the finite longitudinal extension of the source. The 
expansion flow velocity u^(x) is parametrised as 

u^(x)= (cosh rji cosh 7] t (r), sinh?7t(r) e r , sinh^ cosh rj t {r)) , rji—r], r)t(r)=r)f(r/R), (30) 
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with a boost-invariant longitudinal flow rapidity and a linear transverse flow rapidity 
profile. t]f scales the strength of the transverse flow. The scalar product in the exponent 
of the Boltzmann factor can then be written as 



x 



K ■ u(x) = M_|_ cosh(?7 — Y) cosh^(r) — K±— sinhr/^r 



(31) 



Please note that for non-zero transverse momentum K±, a finite transverse flow breaks 
the azimuthal symmetry of the emission function via the second term in (j3~ip . For rjf = 
the source has no explicit K±- dependence, and M± is the only relevant scale. As will be 
discussed in Sec. ^]3| this gives rise to perfect Mx-scaling of the YKP radius parameters 
in the absence of transverse flow, which is again broken for non-zero transverse flow [ [26f] . 

For the numerical calculations below we have selected one fixed set of source parameters: 
R = 3 fm, t = 3 fm/c, At = 1 fm/c, Arj = 1.2, T = 140 MeV. 




800 1000 

MJMeV) 



Fig. 2. Same as Fig.l, but now evaluated in the LCMS. Please note 
the change of sign and magnitude of the cross-term. (Figure taken from 



Ref. [ B7 



6. MOMENTUM DEPENDENCE OF HBT PARAMETERS 

6.1. Standard Cartesian fit 

In Fig. 1 I show the HBT radius parameters from the standard Cartesian fit (|14D for 
pion pairs with cm. rapidity Y = 1.5 where the fit of the correlator is done in the CMS 
[ 27j. The different thick curves correspond to different strengths rjf of the transverse 
flow. Without transverse flow R s is Mj_-independent because the source ( |2"§| ) has no 
transverse temperature gradients. As transverse flow increases, R s develops an increasing 
dependence on M±. It can be approximated by an inverse power law, with the power 
increasing monotonously with rjf [[12|, [14]]. Ri features a very strong M^-dependence even 
without transverse flow, due to the strong longitudinal expansion of the source. It can 
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also be described by an inverse power law, with a larger power 
with the approximate JT/M±-scelmg law suggested in [ 



0.55, in rough agreement 
see, however, [ O, EZTH for a 



more quantitative discussion). The increase of R at small M± is due to the contribution 
flT|) from the effect ive lifetime. As seen in Fig. 4 below, in the YK frame (source rest 
frame) the latter is of order 2.5 fm/c at small Mj_; Fig. lb shows that its effect on R 
compared to R s in the CMS is much smaller (and thus more difficult to measure). Fig. Id 
shows that the cross-term is small in the CMS but non-zero. It vanishes at K± = by 
symmetry and also becomes small again at large K±. 

The thin lines in Fig. 1 show for comparison approximate results for the HBT radii 
calculated from the approximate analytical results given in Ref. [ || which were derived 
by evaluating Eqs. ( |T5] - |I8|) by saddle point integration. It is clear that this method fails 
here (see Ref. [ |12| for a quantitative discussion of this approximation), and that the 
analytical expressions should not be used for a quantitative analysis of HBT data. 
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Fig. 3. (a) The Yano-Koonin rapidity for pion pairs, as a function of 
the pair cm. rapidity Y, for various values of K± and two values for 
the transverse flow rjf. (b) The same, but plotted against K± for various 
values of Y and rjf. (Figure taken from Ref. [ [D|.) 

Fig. 2 shows the same situation as Fig. 1, but now all HBT radii are evaluated in 
the LCMS (longitudinally comoving system [ which moves with the pair rapidity 
Y = 1.5 relative to the CMS. A comparison with Fig. 1 shows the strong reference frame 
dependence of the standard HBT radii. In particular, the cross-term changes sign and is 
now much larger. The analytical approximations from Ref. [ [5| work much better in the 
LCMS [ |U , but for R Q and R\ they are still not accurate enough (in particular in view of 
the delicate nature of the lifetime effects on R D ). 

6.2. The Yano-Koonin velocity 

Fig. 3 shows (for pion pairs) the dependence of the YK velocity on the pair momentum 
K. In Fig. 3a we show the YK rapidity Y YK = ^ln[(l +v)/(l — v)] as a function of the 
pair rapidity Y (both relative to the CMS) for different values of K±, in Fig. 3b the same 
quantity as a function of K± for different Y. Solid lines are without transverse flow, dashed 
lines are for rjf = 0.6. For large K± pairs, the YK rest frame approaches the LCMS (which 
moves with the pair rapidity Y); in this limit all pairs are thus emitted from a small region 
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in the source which moves with the same longitudinal velocity as the pair. For small K± 
the YK frame is considerably slower than the LCMS; this is due to the thermal smearing 



of the particle velocities in our source around the local fluid velocity u^{x) [ [14|. The 
linear relationship between the rapidity Y YK of the Yano-Koonin frame and the pion pair 



rapidity Y is a direct reflection of the boost-invariant longitudinal expansion flow [ [HJ. 
For a non-expanding source Y YK would be independent of Y. Additional transverse flow 
is seen to have nearly no effect. The dependence of the YK velocity on the pair rapidity 
thus measures directly the longitudinal expansion of the source and cleanly separates it 



from its transverse dynamics. A detailed discussion of these features is given in Ref. [ I i 



6.3. Mi-scaling of YKP radii and transverse flow 

In the absence of transverse flow, a thermal source like 



depends on the particle 

2 



rest mass and on the transverse momentum K± only through the combination Mj_ — 
m 2 + K\ (see Eq. (0)). Furthermore, the source is then azimuthally and x — > —x 
reflection symmetric. Hence (xt), (xz), and (x 2 — y 2 ) all vanish and the approximations 
become exact. As a result, all three YKP radii ([22|)- (|24]) are only functions 
as well as of Y, of course), i.e. they do not depend explicitly on the particle 



in Eqs. ([23 
of M]_, too 
rest mass. 
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Fig. 4. The YKP radii R±, Ru , and Ro (from top to bottom) for vanishing 
transverse flow (left column) and for r]f = 0.6 (right column), as functions 
of M±_ for pairs at Y cm = 0. Solid (dashed) lines are for pions (kaons). 
The breaking of the Mj_-scaling by transverse flow is obvious in the right 
column. Also, as shown in the lower right panel, for nonzero transverse 



flow Rq does not agree exactly with the effective source lifetime y(t 2 ). 



(Figure taken from Ref. [ |i4[. 



This is seen in the left column of Fig. 4 where the three YKP radii are plotted for Y cm = 
pion and kaon pairs as functions of Mj_; they agree perfectly. The transverse radius 
here shows no M^-dependence due to the absence of transverse temperature gradients, 
but even with temperature gradients it would only depend on M±. (Of course, this 
discussion neglects resonance decays which will be studied in Sec. |7|.) Note that M±_- 
scaling in the absence of transverse flow applies only to the YKP radius parameters: 
since the expressions QrB|)-(|TBD involve nonvanishing variances with f3±- or /^-prefactors 
(which depend explicitly on the rest mass), the HBT radii from the standard Cartesian 
fit do not exhibit Mj_-scaling. 

For non-zero transverse flow rjf ^ this Mj_-scaling is broken by two effects: first, the 
second term in ( [31]) destroys the Mi-scaling of the emission function itself, and second 



the /3-dependent correction terms in (E3,E3) are now non-zero because the same term 



also breaks, for K± ^ 0, the x — > — x and x —> y symmetries. The magnitude of the 
associated scale breaking due to the pion-kaon mass difference is seen in the right column 
of Fig. 4 for r\f = 0.6. The effects are small and require very accurate experiments for 
their detection. However, the sign of the effect is opposite for R\\ and for R±, Rq which 
may help to distinguish flow-induced effects from resonance decay contributions. 

Since for Y cm = the YK and CMS frames coincide, (3i = in the YK frame and 
the approximation in (1231) remains exact even for non-zero transverse flow. The same is 



not true for the approximation in fl24), and therefore we show in the lower right panel of 



Fig. 4 also the effective source lifetime y (t 2 ) for comparison. The apparently rather large 
discrepancies between the YKP parameter R and the effective source lifetime is due to 
a rather extreme choice of parameters: a large flow transverse flow and a small intrinsic 



source lifetime of Ar = 1 fm/c in (p29|) . Since J (t 2 ) approaches Ar in the limit of large 
Mj_ while the dominant [ [L4| correction term (x 2 — y 2 ) does not depend on Ar, the YKP 
parameter Rq will track the effective source lifetime more accurately for larger values of 
Ar (and for smaller values of rjf). 



Why do y (t 2 ) and Rq increase at small M±? Due to the rapid longitudinal expansion, 
the longitudinal region of homogeneity Ru is a decreasing function M±. Since for differ- 
ent pair momenta Rq measures the source lifetime in different YK reference frames, the 
freeze-out "hypersurface" will in general appear to have different shapes for pairs with 
different momenta. Only in our model, where freeze-out occurs at fixed proper time ro 
(up to a Gaussian smearing with width Ar), is it frame-independent. It is thus generally 
unavoidable (and here, of course, true in any frame) that freeze-out at different points z in 
the source will occur at different times t in the YK frame. Since a z-region of size Rn con- 
tributes to the correlation function, Rn determines how large a domain of this freeze-out 
surface (and thus how large an interval of freeze-out times in the YK frame) is sampled 
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by the correlator. This interval of freeze-out times combines with the intrinsic Gaussian 
width At to yield the total effective duration of particle emission. It will be largest at 
small pair momenta where the homogeneity region R» is biggest, and will reduce to just 
the variance of the Gaussian proper time distribution at large pair momenta where the 
longitudinal (and transverse) homogeneity regions shrink to zero. The rise of At(K) at 
small K is thus generic. 



7. RESONANCE DECAYS 



The proportionality of the Mi-dependence of R± to the transverse flow 77/ and the 
particular pattern of M± scale-breaking by the latter open an avenue for the quantita- 
tive extraction of transverse flow from HBT data [ ^6] . This requires, however, that the 
M^-dependence of R± is not affected by resonance decays. Since they contribute more 
to pions than to kaons they may also affect the M^-scaling arguments. The work by the 
Marburg group [ 10] on resonance decay effects on HBT in the context of hydrodynamical 



simulations indicates, within the standard Cartesian framework and without accounting 
for the cross-term, a possible additional Mj_-dependence of the transverse radius. How- 
ever, a systematic analysis of resonance contributions to HBT as a function of various 
characteristic source parameters is only now becoming available 
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Fig. 5. The influence of resonance decays on the Mi-dependence of R s 
(a,b) and R Q (c,d) for Y cm = pion pairs. a,c: no transverse flow; b,d: 
transverse flow rapidity rjf = 0.3. The Gaussian transverse radius is here 
R = 5 fm, and T = 150 MeV. (Figure taken from Ref. [EH.) 



In Fig. 5 I show some results from Ref. [28] for the same emission function ( p9|) . The 
only change for resonances is an additional spin degeneracy factor and the different rest 
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mass. The complete spectrum of relevant resonances is included, and in the decays the 
2- and 3-body decay kinematics is fully taken into account. The HBT radii are extracted 
from a Gaussian fit to the numerically calculated correlation function. A detailed technical 
discussion is given in Ref . [ |28| . 

Fig. 5 shows that the effects of the short-lived resonances with lifetimes of order 1 
fm/c on R s are essentially negligible, both at vanishing and at nonzero transverse flow. 
Only the uj with its intermediate lifetime of 20 fm/c affects R s , but only for vanishing 
transverse flow. There it induces a weak Mj_-dependence at small M±_ even in the absence 
of transverse flow; at M± > 500 MeV the contribution of the uj dies out, and R s again 
becomes M^-independent (which would not be the case if it were affected by flow). At 
T]f = 0.3 and 0.6 [ Q not even the uj generates any additional M^-dependence! - R Q 
shows some effects from the additional lifetime of the resonances, in particular from the 
long-lived uj. Resonances with much longer lifetimes than the uj (in particular all weak 
decays) have no effect on the radii, because their contribution to the correlator is only at 
very small values of q which cannot be resolved experimentally. They lead to a reduced 



"incoherence parameter" A [ [10], 15]. Since for increasing M± the resonance contributions 



decrease, the A-parameter increases with Mj_, approaching 1 as M± — > oo [ rDj, [TJ|]. A 
detailed study will follow 



The weak effect of resonances on R s = R± seems surprising: due to their non-zero 
lifetime they should be able to propagate outside the original source before decay and 



form a pion "halo" [[TCJ, |15fl . This effect is, however, much weaker than naively expected: 
most of the resonances are not very fast, and the halo thickness is thus only a fraction 
of the resonance lifetime. At finite transverse flow an additional effect comes into play: 
it turns out that then the effective size of the emission function for directly emitted 
resonances is smaller than that for direct pions [|28|]! At i]f=0.3 and 0.6 this even slightly 
overcompensates the halo effect, and altogether the resonances change neither the size 
nor the Mi-dependence of R s . 



8. CONCLUSIONS 

The model-independent expressions of Sees. |] and || for the HBT width parameters in 
terms of second order variances of the emission function provide the basis of a detailed 
physical interpretation of the measured HBT radii. They show that the HBT radius 
parameters do not necessarily measure the full geometric extension of the source, but 
regions of homogeneity in the effective emission function for particles with certain fixed 
momenta. For expanding systems these are usually smaller than the naive geometric 
source size and decreasing functions of the pair momentum. For systems with finite 
lifetime the HBT parameters usually mix the spatial and temporal structure of the source, 
and their unfolding requires model studies. 

With the new YKP parametrization we have found a method which, for systems with 
dominant longitudinal expansion, cleanly factorises the longitudinal and transverse spatial 
from the temporal homogeneity length. The effective source lifetime is directly fitted by 
the parameter R ; it is generically a function of the pair momentum and largest for pairs 
which are slow in the CMS. Another fit parameter, the YK velocity, measures directly the 
longitudinal velocity of the emitting fluid element, and its dependence on the pair rapidity 
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allows for a direct determination of the longitudinal expansion of the source. Without 
transverse expansion, the YKP radius parameters show exact Mj_-scaling. The breaking 
of this scaling and the M^- dependence of the transverse radius parameter R± allow for a 
determination of the transverse expansion velocity of the source. Resonance decays were 
shown to mostly affect the lifetime parameter and leave the M^-dependence of R± nearly 
unchanged. They thus do not endanger the extraction of the transverse flow via HBT. 

With this new and detailed understanding of the method, I believe that HBT inter- 
ferometry has a begun a new and vigorous life as a powerful tool for reconstructing the 
geometric and dynamic space-time characteristics of the collision zone from the measured 
momentum spectra. 
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